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Abstract. A new global basis of B-splines is denned in the space 
of generalized quadratic splines (GQS) generated by Merrien subdi- 
vision algorithm. Then, refinement equations for these B-splines and 
the associated corner-cutting algorithm are given. Afterwards, several 
applications are presented. First a global construction of monotonic 
and/or convex generalized splines interpolating monotonic and/or con- 
vex data. Second, convergence of sequences of control polygons to 
the graph of a GQS. Finally, a Lagrange interpolant and a quasi- 
interpolant which are exact on the space of affine polynomials and 
whose infinite norms are uniformly bounded independently of the par- 
tition. 



§1. Introduction and notations 

This paper is a continuation of [7] and [10] where a one-parameter family of 
C 1 Hermite interpolants was defined by the Merrien subdivision algorithm 
([6], abbr. MSA). This family provides a general solution to arbitrary 
monotone and convex Hermite interpolation problems with data at the 
ends of a single interval. Here we extend and complete these results by 
considering Hermite interpolation on a partition of a given interval by C 1 
functions whose restrictions to subintervals are generated by MSA. Let 
X := {a = xq, x\, . . . , x n = b} be such a partition of / = [a, b] with 
Ii = [xi-i, Xj\ and hi = Xi — x^-\ for 1 < i < n. Let (3 = (fa, . . . (3 n ) be a 
sequence of parameters fa G [— 1,0[. 

On each subinterval Ii, we consider the 4-dimensional space V(fa) gener- 
ated by Merrien subdivision algorithm and depending on the specific pair 
of parameters (di,fa), where aij = wxzjf) e [ — ^ , 0[. In that case, it was 
proved that the MSA is C 1 -convergent. Moreover, for / G V(fa), f is 
Holder and satisfies the following inequality for some constant C > and 
for ji = —log2(l + \fa) (see [6], proposition 2): 



\f(x)-f(y)\<C\x-yP\ (x,y)elixli. 
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We denote by GSzil, X, (3) or simply GS2HS) the subspace of all functions 
g G C 1 (J) whose restriction to Ii is in V(f3i) for 1 < i < n. The dimension 
of this space is 2n + 2 since g is uniquely determined by Hermite data 
Hi = g(xi) and pi = g'(xi) for < % < n. The elements of GS2(f3) are 
called generalized quadratic splines (abbr. GQS) because their properties 
are similar to those of classical quadratic splines, which correspond to the 
choice f3i = — 1, 1 < % < n. 

Here is an outline of the paper. In section 2, we define a new global basis of 
B-splines for the space GS2{(3) of GQS. In section 3, we use these B-splines 
to express in a new global form the results on monotonicity and convexity 
preserving properties of GQS, already given in [7] [10] for functions defined 
on one subinterval. In section 4, we also give a global version of the 
algorithms constructing monotonic or convex interpolants given in the 
same papers. When compared to other C 1 Hermite interpolants of the 
literature, the advantage of our method is its simplicity both in the local 
construction of the GQS and in the adaptivity of the method to arbitrary 
sets of data. For example, in [1][2], the method is adaptive, however it 
may need polynomials of arbitrary high degree, which can lead to rather 
complicated expansions and calculations. 

In section 5, we give the refinement equation for coarse B-splines in terms 
of fine B-splines in the space of GQS defined on a refinement of the initial 
partition. From this result we deduce a corner-cutting algorithm (see [9] 
for definitions and properties) which is the geometric form of the algorithm 
expressing the new coefficients of a GQS g in the B-spline basis of the fine 
space in terms of its old coefficients in the B-spline basis of the coarse 
space. We also prove the convergence to the graph of g of the sequence of 
control polygons associated with successive steps of this algorithm. 
Finally, in section 6, we study two approximation operators: a Lagrange 
interpolant and a quasi-interpolant which are both exact on the space 
Pi of affine polynomials and whose uniform norm are uniformly bounded 
independently of the given partition on i". This extends a previous result 
given for ordinary quadratic splines by Kammerer, Reddien and Varga [4] 
and also by Marsden [5] . We postpone numerical applications to a further 
paper which should contain variants and refined versions of the general 
algorithms presented in Section 4 of the present paper. 

For the sake of clarity, we now recall the basic equations of the MSA giving 
the values at the midpoint m = |(a + b) of [a, b] of a function / and its 
first derivative /' from the four values: 

{/(a), /'(a); /(&),/'(&)} 

at the endpoints of the interval (see [5]). The construction starts with 
[a, b] = Ii = [xi-i, Xi] and gives the values of / and /' at the dyadic points 
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of Ii. Let h = b- a and 0* = ^^ri = -2«; e]0, |], then 



/("*) = ^ ((/(«) + /(&)) " W(b) - f'(a))) 

f'(m)- 1 f iJj 2fl f'(b) + f'(a) 
f (m) " 1^207 V & 1 2 

In each subinterval lj, let us define the two points 

£i = Xi-6ihi, rji = Xi + 9 i+ ihi+i, 

with = an d r]n = x n . Then each element ^ G V(/%) can be expressed 
as 

gi = ai-xbo + + Cib 3 + 

in the local B-spline basis {bo, b\, &2> ^3} of V(0i) defined in [10]. By 
definition, the quadruplet [eij_i, Q, a$] is the list of B-coefEcients of 
<7j on the subinterval 7^. The iocal control polygon (abbr. LCP) of gi has 
the four following local control vertices 

Oi_i = (xi_i,ai_i), = (7/i_i,di_i), Cj = (Cj,Cj), a i = (^,a i ). 

The ordinates of these points are the B-coefficients of gi and they can be 
expressed in function of the four Hermite data (yi-i,Pi-i',yi,Pi) at the 
two end-points of If. 

Vi = 9i(xi) = 9i+i(xi), Pi = g'i(xi) = g' i+1 (xi). 

They are given, for < % < n — 1, by 

a>i = Vi, Ci = at - OihiPi, di = a t + 9 l+1 h i+1 pi, a i+1 = y i+1 

(with the convention ho = h n +\ = 0). 

Using the MSA and the properties ol\ = —\0i, A = 20—1 ■> 1 ~~ A = 
1 _ 26) , and £j — 77^—1 = (1 — 29f)hi, we obtain respectively at the midpoint 

Xi — l + Xi 

rxii = of L: 



1 1 

9i{™*) = g (( a *-l + a _ ^i^i(P< - Pi-l)) = -^{di-1 + Ci), 

// \ 1 fcti-di-! nQ Pi-i+Pi\ Ci-di-x 

= (— — M '^^J = 173^-7. 

which proves that the tangent to the curve at rrii is the segment rfj_iCj 
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§2. Generalized quadratic B-splines 

In this section, we define a global B-spline basis of the space GS 2 (f3) of 
generalized quadratic splines. 

Let Ui — — — for 1 < i < n — 1. 

Villi + tfj+l/lj+l 

Definitions 

(i) For 1 < i < n — 1, let B 2i be the function whose support is the segment 
[xi-i,mi + i] and whose lists of B-coefficients in the subintervals Ii and 
Ii+i are respectively 

[0,0,1,^], 0,0,0]. 

(ii) Similarly, let B 2 i + i be the function whose support is the segment 
[rrii, Xi+i], and whose lists of B-coefficients in the subintervals U and U + i 
are respectively 

[0,0,0,1 -a;*], [1-^,1,0,0]. 

(iii) Moreover, there are 4 special B-splines {B , Bi, B 2 n, B 2n +i} at the 
end points of J, defined respectively by 

supp(-Bo)=[xo, mo], its list of B-coefficients in I\ being [1,0,0,0]. 
supp(-Bi)=[xo, xi], its list of B-coefficients in I\ being [0,1,0,0]. 
sup-p(B 2n ) = [x n -i, x n ], its list of B-coefficients in I n being [0,0, 1,0]. 
supp(S2 n +i)=[m n , x n ], its list of B-coefficients in I n being [0,0,0,1]. 

It is easy to verify that 

B 'oM = V = ^T~' B 2n(Xn) 



B[(x ) = -B' (x ), B' 2n+l (x n ) = -B 2n (x n ), 
and for 1 < i < n — 1 

BU*i) = —V = -(TIT- = 

Theorem 1. The generalized quadratic B-splines {B k ,0 < k < 2n + 1} 
form a basis of the space GS 2 (J3). Moreover, they form a partition of 
unity, or a blending system, in this space. 

Proof: First, let us prove that the B-splines belong to the space GS 2 {J3), 
i.e. that they are C 1 continuous at the points {x\,x 2 , . . . ,x n -\}. By 
construction, they are already C 1 in each subinterval Ii. In addition, 
their derivatives satisfy the above relations, so they are continuous at the 
interior points of X. 
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The B-splines are linearly independent: assume that g = Ylk^) 1 IkBk = 0. 
Then, the B-coefficients of the restriction gi = XV=^2i-2 IkBk of g to Ii 
are respectively 

O-i-l = ^t-llH-l + (1 - ^i-l)72i-l) di-l = 72i-l) 

Q = 72i, a* = ua2i + (1 - a>i)72i+i- 

Therefore, since a;j_i and 1 — cu^ are non zero, we get successively 724-1 = 
72i = 0, and 72^-2 = 72i+i = 0. 

The B-splines generate the space GS 2 ((3): it suffices to express the co- 
efficients 7fc in function of the B-coefficients. The above equations give 
immediately 

72i = Ci, 724+1 = di. 

Finally, let us prove that J^-Bfc = 1: it suffices to prove it on each 
subinterval 1^. As local B-spline bases are blending systems, one has 
to show that the sum of local B-coefficients of global B-splines is equal 
to 1. This property is easily deduced from the lists of B-coefficients of 
B 2 i-2, B 2 i-i , B 2i , B 2 i+\ on the interval Ii which are respectively equal to 

[u>i_i,0,0,0], [1 1,0,0], [0,0,1, u^], [0,0,0,0;*] 

□ 

§3. Global and local control polygons. Monotonicity and convexity 

It has been proved in [6] that the spaces V(/%) contain the space Pi 
of affine polynomials . As the list of local B-coefficients on Ii of the 
function e\{x) = x is [a^-i, ?7i_i, Xj], the results of section 2 show that 
the representation of e\ in the basis of generalized B-splines is given by 

n+l 

ei = ^2(^iB 2i + r]iB 2i+1 ). 

i=0 

Therefore, we can define the global spline control polygon (abbr. SCP) of 
g = Efc=o l^Bk with vertices 

72i = (&»72i), 72i+l = (?/i,72i+l) 

from which we easily deduce the vertices of the local control polygon (LCP) 

Cli-l = <jJi-i^2i-2 + (1 - Wi_i)72i-i, Ci-i = 72i-2, 
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di = 72i+l, CLi = UJij2i + (1 - UJi)^2i+l- 



Theorem 2. A function g e GS 2 ((3) is monotonic (resp. convex) if and 
only if its SCP is monotonic (resp. convex), with the same sense of varia- 
tion. 

Proof: ^From theorem 6 of [10], we know that g is (e.g.) increasing on Jj 
if and only if its LCP is increasing, i.e. iff 

CLi-l < di-i < Ci < di. 

Since we have 

di-i - a,i_i = cui_i(72i-i - 724-2), Ci - di-i = 72i - 724-1, 

O-i - Ci = (1 - ^i)(72i+l - 72i), 

we see that the above inequalities are satisfied iff 724-1 < 724-1 < 724-1, 
i.e. iff the global SCP is increasing. 

Similarly, from theorem 8 of [10], we know that g is convex on Ii if and 
only if 

di-i — o>i-\ Ci — di-i di — Ci 

Vi-l - x i-l ~ Ci - Vi-l ~ Xi~£,i 

Using the equalities 

Vi—i si— 1 = , Si - Vi-l = I 1 - 20i)fti, Vi-i. 



UJi-i 1 — u 



% 



we obtain successively the following identities 

di-i — di-i _ Ui-i (724-1 — 724 - 2) _ 724-1 — 724-2 
Vi-l - x i-l &ihi Vi-l ~ Ci-1 

Cj - dj-1 _ J2i - 724-1 _ 724 - 724-1 

£4-774-1 (l-20i)hi &-Vi-i 
o-i - Ci (1 - a;,) (724+1 - 72*) 724+1 - 724 



x i 4*4 @ihi Tji ^i 

which show that the convexity of the global SCP is equivalent to that of 
the LCP in each subinterval, i.e. to the convexity of g. □ 



§4. Construction of monotone or/and convex interpolants 
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Algorithm 1: monotone interpolant 

This algorithm describes the construction of an increasing interpolant 
g £ GS2(I,X,{3) to arbitrary increasing Hermite data: 

(Vi,Pi), 0<i<n, 
assumed to satisfy the properties: 

A Vi > 0, Pi > 0. 

(There is a similar algorithm for decreasing data). The main point consists 
in choosing the sequence (3 in function of the data. ^From theorem 2, we 
know that g = J2 k jkBk is increasing if and only if the sequence (7^) of 
its S-coefhcients is increasing. On one hand, we already know that 

72i+i - J2i = di- Ci = (6 i+1 h i+ i + 9ihi)pi > 0. 

On the other hand, using the notations: 



we must have 

72; - 72i-i = Ci - di-t = hiiji - 26iHi) > 0, 
which implies the following condition on the parameter 9i = \(3i(3i — 1 : 



< 6i < K = 




Note that this condition is equivalent to g'irrii) > 0. 

There appear two cases and we obtain the following 

Theorem 3. (i) if fa < 2T i} then we can choose 6i = \: in that case, 
the local interpolant in the subinterval Ii is an ordinary quadratic spline 



case, the local interpolant in the subinterval U is a generalized quadratic 
spline associated with the parameter fa = r 2 J . 



j. 



-» Hi = o(Pi-l +Pi) 



(Pi 




< \: in that 



8 



ps 



Algorithm 2: convex interpolant 

This algorithm describes the construction of a convex interpolant 
g £ GS2(I,X,/3) to arbitrary convex Hermite data: 

(y»Pi), 0<i<n, 
assumed to satisfy the properties 

< n < pi 

on each subinterval U. (There is a similar algorithm for concave data). 
Using the results given in the proof of theorem 2, we get the convexity 
conditions on Ii for the parameter 6% = \ — 1 : 

Ti - 26im 

Pi -' - jr^y - Pi 

which can also be written 

(1 - 6i)pi-i + QiPi < n < (1 - 9i)pi + 9iPi-i 

As for algorithm 1, there appear two cases: 
Theorem 4. 

(i) if Ti satisfies the two inequalities 

j(3pi-i + Pi) < n < \{pi-i + 3pi), then we chose 9i = \: in that 
case, the local interpolant is an ordinary quadratic spline. 

(ii) else the local interpolant is a generalized quadratic spline (0$ < |). 

a) either pi-\ < Ti < |(3pi_i +Pi): then we have to choose 9i < 9~i = 

Tj-Pi-I 
Pi-Pi-! ' 

b) or \{pi-i + 3pi) < Ti < pi\ then we have to choose 9i < 9~i = 

Vi-Tj 
Pi-Pi-! ' 

In both cases, the GQS belongs to V with fa = {^e- • 

Algorithm 3: monotone and convex interpolant 

For sake of simplicity, we assume that the data are increasing and convex 

< p^ <Ti <Pi 

on each subinterval Ii. (There are similar algorithms for the three other 
cases). By putting together conditions of theorems 3 and 4, we get the 
same algorithm as in the preceding case. 
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Theorem 5. (i) if |(3pj_i + Pi) < Ti < ^(pi-i + %Pi), then we can 
choose 9i = \. The local interpolant is an ordinary quadratic spline. 

(ii) else the local interpolant is a generalized quadratic spline. 



a) either Ti < j(3pi-i +Pi), then we chose 9i < 9i 



Ti—Pi 



Pi-Pi-i 



b) or n > \{pi-i + 3pi), then we chose 9i < 9 t = Pi Ti . 

In both cases, the GQS belongs to V($), with fa = ■ 
Remark: these are only rough algorithms: in practice, one has to smooth 
a little bit the above conditions and also to treat the cases when there 
appear equalities in the conditions. This will be done in a further more 
complete paper illustrated with numerical examples. 

algorithm ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



§5. Refinement equations and global corner-cutting algorithm 

In this section, we consider the subpartition X = X U {m 8 , 1 < i < n} 
dividing / into 2n subintervals. The space GS 2 (I,X, 0) is a subspace of 
dimension 2n + 2 of the new space GS^I, X, /3) of dimension 4n + 2. Here 
j3 denotes the sequence of parameters deduced from (3 by taking twice the 
same parameter 0, once for the left subinterval 1[ = [xi-i,rrii\ and once 
for the right subinterval I'/ = [m^a^] of If. 



The finer B-splines of GS 2 (I, X,(3) are denoted 



{B h < / < 4n + 1}. 



Thanks to the local corner-cutting algorithm (abbr. CCA, see [8], section 
), one can compute the local B-coefficients in the subintervals I[ and I'/ 
of the coarser B-splines {Bk, 0</c<n + l}in function of their previous 
B-coefficients in Jj, and we get the following refinement equations. 

Let us recall that Ui = 1+1 1+1 , for 1 < i < n — 1, uo = l — uj n = 1. 

9ihi + 9i + ihi + i 

Theorem 6. (i) For all 1 < i < n— 1, one has the two following refinement 
equations: 

11- 11- 1 1 

B 2 i = (- + -^Pi)Bu-2 + (- - -^(3i)Bu-i + -(1 + Ui)B u + -UiB 4i+1 , 

B 2 i+1 = UJ i)Bu+ -^(2— UJi)B 4: i + i + (- — -{3i)B4i + 2 + (- + -fli)B4:i+3- 



10 ps 
(ii) For the B-splines at the endpoints, one has respectively: 

B = B + ISx, B x = l -B x + (i - i/^OBa + (1 + 1^)53, 

B 2n = 2"^ 4n "'"( 2 ~~ ~^Pn)B4n-l + ( ^ + ~~^fln)B4n-2, i?2n+l = ^ B^n + B^n+l ■ 

Proof: We only give the proof for B 2 i, the others being similar. Let 

B2i = ^u-2Bu-2 + \i-u-\Bu-\ + fluBu + H4i+lBu+l- 

By application of the CCA, starting from the B-coefficients [0,0, l,u)i] of 
B 2i on the interval 7^, we deduce the B-coefficients of B 2i , respectively on 
the subintervals I[ and I'/: 

[0,0, \ + \^\\ and + 

Similarly, from the B-coefficients [u>i, 0, 0, 0] of B 2 i on Ij+i, we deduce the 
B-coefficients of B 2i on the subintervals and 

[^,^,0,0] and [0,0,0,0] 

On the other hand, the B-coefficients of the finer B-splines are respectively 

1) on the subintervals I[ and I" 

for B 4i _ 2 : [0,0,1,^], [^,0,0,0] 

for B u _ x : [0,0,0,^], [^,1,0,0] 

2) on the subintervals l'{ and 

for B u : [0,0,1,^], K 0,0,0] 
for B 4l+1 : [0,0,0,1 -Wi], [1-^,1,0,0] 

Therefore, the B-coefficients of B 2 i as linear combination of the four finer 
B-splines on the three subintervals I[, I" and are respectively equal 
to 



[0, 0, flU-2, 2 (A*4i— 2+A*4i— l)] j (A*4i-2+A t 4i-l), /Mi-l, ^4i, Wi^4i + (1 — Wi)/X4i+l] 
[c^/i 4 ; + (1 - 0;i)/i4i+l, A*4i+1,0,0]. 

By identifying these B-coefficients with those of £?2i, one obtains: 
11. 1 1. 1,, v 1 

^4i-2 = ^ + A*4i-1 = - - ^Pi, ^4i = "(I /X4 i+ i = -UJi 

□ 

An immediate consequence of the previous theorem is the following global 
corner-cutting algorithm: 
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Theorem 7. Given the two expansions of S G GS^i", X, (3) with respect 
to the coarse and fine B-splines bases: 

2n+l 4n+l 
k=0 1=0 

then, the new B-coefEcients have the following expressions in terms of the 
former B-coefEcients: 

= (\ + |A)72i-l + (\~ 
hi = 7j(l + ^i)72i + ^(1 - UJi)j2i+l, 
= ^il2i + 7^(2 - CUi)7 2i+ l. 

Proof: The proof simply consists in comparing the coefficients in the two 
expressions after substituting in the first expression the coarser B-splines 
Bk by their expansions as linear combinations of the finer B-splines B\, 
given in theorem 3. □ 

Theorem 8. The sequence of SCPs associated with a given GQS, ob- 
tained by successive applications of the global CCA, converges uniformly 
to the GQS. 

Proof: Let S(x) = $^fc="o IkBk be the equation of a GQS and let Po = 
Sfc^o 1 lk<ftk be the equation of its initial SCP. We denote by 4>2i the hat 
function with support [77^-1, 7/j] satisfying 02i(&) = 1) an d by <p2i+i the 
hat function with support satisfying 4>2i{j]i) = 1- 

For x G [r]i-i,€i], (resp. x G [£;,??;]), we have 

\S(x) - P (x)\ < max{\S(x) - 72i _i|, \S(x) - 72l |} 

(resp.\S(x) - Pq(x)\ < max{\S(x) - 72* |, \S(x) - 72^+1 1}). 
Moreover, for x G [f7i-i,£i], we have 

2i+l 

|5(x)-72i|< l7fc-72i|-Bfc(x) < max{|7fc-72i|,2i-2 < k < 2i + l} 

k=2i-2 

As we observe that 

\lk - 72*1 < 2max{|7 fc+ i - 7 fc |,2i - 2 < k < 2i + 1}, 
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we are led to define 

A = max{ |7fc+i - 7fc|,0 < k < 2n}, 
and we finally obtain 

\\S — -Pol |oo < 2Ao- 

Now, we do the same for the next SCP Pi deduced from P by one appli- 
cation of the global CCA. Here, we define 

Ai = max{\5i +1 - 5j|,0 < I < 4n + 1}. 

We have successively the following majorations 

|<J 4i _i - 6u- 2 \ = I 2 A (72< - 72i-i)| < ^l72i - 72i-i| 

1^47+1 - &ti| = - |72i+l _ 72i| 
|<5 4l - 5 4l -i| = |^(1 - Wi)(72i+1 - l2t) + (\ + ^A)(72t - 72i-l))| 

< ^ max{|7 2i+ i - 72i|, (72* - 72i-i|} 
Therefore we obtain 

Ai < ^A . 

Denoting by A m the maximum distance between two consecutive vertices 
of the SCP P m obtained after m applications of the global CCA, we get 

1 A 

A m < — A , 

which proves the uniform convergence to the GQS g of the sequence (P n ) 
of its SCPs. □ 



§6. Quasi-interpolant and Lagrange interpolant 

In this section, we define two kinds of approximation operators: the first is 
a quasi-interpolant with good shape-preserving properties while the second 
is a Lagrange interpolant which is uniformly bounded independently of the 
partition. Let us begin with the quasi-interpolant Q defined by 

n+l 

Qf = ^[/(£i)£ 2i + f(vi)B2 i+ i] 

i=0 
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Theorem 9. Q is exact on Pi, ||Q||oo = 1, and Q preserves the mono- 
tonicity and the convexity of f. 

Proof: Since e = J2™=di B 2i + B 2i +i] and e x = YJl=o\ii B 2i + Vi B 2i+i], 
we obtain immediately the property that Q is exact on Pi. Moreover, 
Qe = e and |Q/|oo < \f\oo imply ||Q||oo = 1- 

If / is increasing (resp. convex), the global B-polygon of Qf is also in- 
creasing (resp. convex), and the result follows by a direct application of 
theorem 2. □ 

Now, let us study a Lagrange interpolation operator L in the space GS2 (P) • 
Let m\ (resp. m") be the midpoint of I[ = [xj_i,mj] (resp. of 1" = 
[rrii,Xi\. The following theorem is an extension to GQS of a result pre- 
viously given by Kammerer et al. [4] and by Marsden [5] for ordinary 
quadratic splines. 

Theorem 10. (i) Given a function f defined on I = [a,b], there exists 
a unique generalized quadratic spline Lf e GS2(I,X, (3) satisfying the 
following interpolation properties: 

Lf(a) = f(a), Lf(b) = fib) and for 1 < i < n, 

Lf(m'i) = /(mj) and Lf{m'() = f{m'(). 

(ii) The B-coefEcients of Lf = Ylt^Q 1 lk B k solutions of the tridiagonal 
system of2n linear equations (1 < i < n): 

a;j_i72i_2 + (3 - c^_i - ^A)72i-i + (1 + ^ Pi)l2i = 4/(m-) 

(1 + - A)72i-i + (2 + Ui- ^i)l2i + (1 - Wi)72i+i = 4/(m") 
with 70 = f(a), 7 2n+ i = f(b). 

Proof: On the interval Ii, we are able to compute the values of Lf = 
Ylt l =2i-2 IkBk at the points and m" by applying the CCA to B-splines. 
Without going into details, it is straightforward to obtain the coefficients 
of the two equations given in the theorem. The tridiagonal matrix of the 
linear system is strictly diagonally dominant since the two corresponding 
inequalities 

3 - u t -t - > iVi-x + 1 + and 2 + u>i - ^fa > 1 + ^ + 1 - ^ 
are respectively equivalent to the following 



2(1 - Ui) - Pi > and 2u>; - Pi > 
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which are obviously satisfied since /5, < and < Ui < 1. □ 

Theorem 11. ||L||oo is uniformly bounded for all partitions of I. More 
specifically, setting j3 = max{fli, 1 < i < n}, we obtain 

4(30-1) 



\L\\oo — 



0(5-30)' 



Proof: For this purpose, we start from the expression of Lf = J2 r =o a i,rW 
in the local B-spline basis of each subinterval J^. Then, by the CCA, we 
compute the values of the b r at the two points m! i and m'( and we get the 
two equations: 

4L/K) = a,, + (2 - \$i)a^ + (1 + ^K, 2 = 4/(mJ) 

4L/(mJ') = (1 + ^A)a M + (2 - ^AK,2 + a l)3 = 4/(m<) 

We have to add to these equations the C 1 continuity conditions of Lf at 
interior points Xi, which can be written: 

«i-i,3 = cii,o and Ai_i(ai_i )3 - ai - 1, 2) = A^a^i - a i)0 ) 

where A l = 2 (^p) = ^- for all 1 < i < n. 
^From these conditions, we deduce the expressions 



Oi-1,3 = «i,0 = 



Ai-i + Ai 



that we substitute in the two equations above. Taking as unknowns 
a,2i-\ = a^i and a 2 i = a^2, we then obtain the following tridiagonal 
system of equations 



A 



i-l 



\-i + A, 



■o-ii-i + 



(1 + -A)a 2i _! + 



Ai-i + Ai 



o 2 i_i + (1 + -0i)a 2i = 4/(mJ) 



Ai + Ai_|_i 2 



&2i + 



Ai+1 a n/ //\ 

-a 2i +i = 4/(m i ) 



Ai + Aj+i 



The tridiagonal matrix is strictly diagonally dominant since the two con- 
ditions 



A,- 



Ai-i + A, 



Ai 



Ai + Ai_|_i 2 2^ Ai + Ai_|_i 



Generalized quadratic B-splines 
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can also be written 

~Pi > 



-2X, 



and — f3 i > 



-2A, 



\-i + \ \+i + \ 

and they are obviously satisfied since (3i < and Aj > for all %. 
Let I 

fl|oo — max{\ak\, <J k <J 2?t, + 1}. ^From the above tridiagonal 
system, we deduce respectively the following inequalities 



Ai-i + Aj 

A ? ; 



low-il <4|/K)|+ 



A, ; - 



Ai-i + A,- 



\a 2 i-2\+(l+^(3i)\a 2 i\, 



Xi + Aj+i 2 ^ 

and we obtain for all i: 

2A ?; 



\a 2i \<A\f{m'l)\+{\+ l -p i )\a 2i _ 1 \+ k+1 



Xi + Aj+i 



I a 2i+l|, 



Ai-i + Aj 



-A 



<4||/|| 00 and 



2A, 



Ai+i + A, 



-A 



a 



oo <4 



Defining A = max{Ai, 1 < % < n} and (3 — max{Pi, 1 < i < n}, then we 
have 4 < A, < A = 2^, hence 



2A 7 - 



Ai_i + A, ^ - A,_i +4 
Therefore, we finally obtain 

4(3^-1) 



A + 4 



= 



1(5 - 3g) 
(3)9 - 1) 



< 



/3(5 - 3/3) 

As the graph of Lf lies in the convex hull of its local control polygons, it 
is not difficult to see that we have also 

4(30-1) 



< 



(3(5-3(3) 

In other words, we obtain the following upper bound which is independent 
of the given partition, but only depends on the sequence (3: 

4(3)9-1) 



oo < 



0(5 spy 



□ 



Remarks: 1. In the case of classical quadratic splines, we have (3 = — 1, 
whence ||£||oo < 2, which is the sharp upper bound given in [4] and [5]. 
This suggests that the upper bound given in theorem 11 may also be sharp. 

2. As for error bounds, we know, from a classical result in approximation 
theory (see [3], theorem ), that for P = Q or L 

\\f - Pf Woo < (1 + ||P||oo)doo(/, GS 2 (p)). 
For / G C 2 (I), we already know (see [3]) that d^f \GS 2 (p)) = 0(h 2 ) , 
with h = max hi. Therefore, we see that both operators have an approxi- 
mation order equal to 2. 
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